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ABSTRACT 

Saturn IB/Centaur and Atlas/Centaur launched electric propulsion vehicles 
employing present or near-future state-of-the-art power supplies (50 to 
100 lh/kW e ) are shown to he attractive systems for performing the 0. 1-AU solar 
probe mission. This result compares favorably with the alternative systems 
composed of the same boosters but with chemical or nuclear final stages - which 
are inadequate for this mission. The performance calculations for the electric 
propulsion systems include the optimization of the launch vehicle burnout 
velocity, the electric stage specific impulse and power level, and the thrust 
orientation for constant specific impulse thrustors. Two sets of results are 
presented. The first set is associated with constant power operation of the 
thrustors as would be the case for nonsolar dependent power supplies such as 
nuclear- electric systems. The second set of data represents solar- cell- 
powered electric systems. In this case the thrustor power is a function of 
the sun-vehicle distance. The power is varied by adjustment of the propellant 
flow rate. Of special interest is the fact that absolute optimal trajectories 
for the solar cell systems belong to an entirely different class than those 
optimals associated with constant power systems. 


*Aerospace Research Engineer 
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initial acceleration of electric stage, ft/sec 2 
power profile function 
32.174 ft /sec 2 

specific impulse of chemical stage, sec 
specific impulse of electric stage, sec 
hardware fraction of chemical stage 

ratio of tankage to propellant mass of electric stage 
hardware mass ' {propellant tanks, engines, structure, etc.) of chemical 
stage, Tb 

initial mass of vehicle in circular Earth orbit, lh 
propellant mass of chemical stage, lh 
final mass, lh 
payload mass, lh 

initial mass of electric stage, lh 

propellant mass of electric stage, lh 

powerplant mass of electric stage, lh 

structure mass of electric stage, lh 

tankage mass of electric stage, lh 

total power delivered to electric thrustors, kW 

heliocentric radius, AU 

initial circular Earth- orhit radius, ft 

radius of sphere of influence, ft 

chemical stage humout velocity, ft/sec 

circular Earth-orhit velocity, ft/sec 
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esc ci^jG velocity (value of V |2 
initial heliocentric velocity, ft/sec 
velocity at r s , ft/sec 
Earth’s heliocentric velocity, ft/sec 

overall specific mass of electric powerplant (including power condition- 
ing and ion engine system) at 1. 0 AU, lb/kW e 
thrust or efficiency 

gravitational constant of the Sun, ft^/sec^ 

terminal configuration angle between Sun-probe line and Sun-Earth line, 
deg 

Subscript 

optimum.' ::, o ' v; 

INTRODUCTION 


Inasmuch as the Sun is unquestionably the very basis of our solar system 
and since it is the predominate factor controlling the environment of the 
solar system, it is quite important that new scientific knowledge be found 
about this "center of attraction. " However, Earth-based solar studies are 
severely hampered by the protective atmospheric and electromagnetic shields 
surrounding the Earth. Of course, seme solar study has been and will continue 
to be carried out by interplanetary probes and orbiting solar observatories. 

But several of the important proposed solar experiments (corona radar sounding, 
ultraviolet and X-ray spatial resolution of the solar disk, determining the 
connection between the rotating solar corona and the solar wind, etc. ) require 
close-in solar probes; that is, probes whose perihelion radii are on the order 
of 0.1 AU. This paper analyzes vehicle systems capable of accomplishing 
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this mission (with emphasis on the solar- cell-powered electric system) rather 
than the solar experiments or the scientific payloads required to perform them. 

The major conclusion regarding all high- thrust propulsion for the 0.1 AU 
solar fly-by probe mission 3 3 is that only very small payload fractions 
(if any at all) are realizable. In particular, both the all chemical : 

Saturn IB/Centaur/kick stage and the Saturn IB/nuclear (l-i stages) launch 
vehicles cannot deliver any payload on a 0. 1 AU perihelion radius trajectoiy. 

It is true that Saturn V boosters with additional high-energy chemical upper 
stages or a nuclear upper stage can forcibly deliver some payload for this 
mission, but the payload fractions are still unattractively small and such 
systems are at least a decade from becoming operational. 

As is usually the case when high-thrust propulsion capabilities are pressed 
close to their limits, interest in electric-propulsion systems is stimulated. 
Past electric-propulsion studies® ' i:0 ® for this mission have, however, assumed 
rather advanced technology - either by choosing quite low values (e.g. , 10 to 
30 lb/k_W e ) of the overall specific powerplant mass a or by assuming a large 
Saturn V/nuclear- electric system. A recent report® by the present author 
presented variational results for this mission that showed that chemical- 
electric hybrid systems offered attractive capability even with present and 
near future power-supply technology (i.e. , Saturn IB/Centaur/electric vehicles 
were assumed, with a = 10-150 lb/kW e ). These results were generated with 
electric power held constant. Other recent studies- 1 - 0,11 have shown that 
lightweight solar- cell power supplies can be manufactured using current 
technology with a's on the order of 50 lb/kW e . 
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Since the power output of a solar cell varies with the distance from the 
sun, the results of Ref. 9 are not applicable to the case of solar- cell 
powered probes. The constant power optimal trajectories of interest reported 
in Ref. 9 would cause the vehicle to spend the major portion of its trip at 
considerable distances outside of the Earth’s orbit. Hence, solar-powered 
vehicles might be expected to suffer severe performance penalties when compared 
with like constant power vehicles. With the intent to clarify this, the present 
paper concerns itself primarily with solar- cell-powered solar flyby probes. 

Also, the great potential of electric propulsion for the 0.1 AU- solar- 
probe mission by using current state-of-the-art technology is demonstrated. 
Naturally, a multitude of systems integration problems should be considered 
for any new type of spacecraft design. A survey of such problems is contained 
in Ref. 12. An in-depth investigation into such design areas as accommodation 
of the changing solar array output to the engines, and deployment and orienta- 
tion of the solar array for a Mars orbiter mission by using solar electric 
propulsion is reported in Ref. 13. 

ANALYSIS 

The analysis presented below is concerned with determining the performance 
of electric-powered solar flyby probes. Rather than repeat analyses already 
reported for chemical and nuclear powered systems, results for these systems 
will be taken from other studies as required. The primary performance 
criterion is the payload ratio. The electric stage is assumed to start either 
in an initial circular 100-nautical-mile Earth orbit or in heliocentric space 
after a chemical boost out of the reference orbit. These two cases are 
referred to as the all electric and chemical- electric hydrid systems, respee- 
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tiveiy. The optimization of the amount of chemical velocity audition "beyond 
the reference point would require solutions to the three-body variational 
problem whenever the chemical boost velocity were less than escape velocity. 
Such solutions are quite difficult to obtain from a computational standpoint; 
therefore, the chemical boost velocities were restricted to values at least 
as large as escape velocity. 

Chemical-electric hybrid system . - The mathematical problem for the hybrid 
system may be stated as follows: given a mission time and specific-power 

powerplant mass, determine the maximum payload ratio achievable with a 
continued firing of the initial launch vehicle and an electric final stage by 
finding the optimum trajectory, chemical stage burnout velocity, electric 
stage specific impulse and initial acceleration. In order to avoid the complex 
calculations involved in the three-body variational problem as discussed here, 
three assumptions were made: (l) the chemical stage is fired until at least 

escape velocity is attained; (2) the gravitational effect of the Sun is 
negligible during the chemical stage of the flight, while the gravitational 
effect of the Earth is negligible during the electric stage of the flight, 
and (3) the time elapsed between chemical stage burnout and sphere of 
influence penetration is negligible. 

The thrust vector control of the heliocentric phase is determined by 
variational principles. This statement means that (l) the Euler- Lagrange 
equations are employed for the determination of the electric thrustor orien- 
tation, (2) the transversality relations are used to optimize the heliocentri c 
travel angle and Earth escape velocity orientation, and (3) the coast phases 
are optimized. 
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The payload equation for the chemical-electric hybrid system just 
described may be written 

= ^ r _ ^ _ ^pp _ . <\ (1) 

Mq M 0 ^ m 0 m 0 m Q m 0 J ' 

where m o /M 0 is the payload ratio of the chemical stage from Earth orbit and 
the sum of terms enclosed in parentheses represents the payload ratio of the 
electric stage. In order to facilitate the discussion concerning the maxi- 
mization of the payload ratio ®l/m o equation (l) will be rewritten in terms 
of the pertinent problem variables (specific impulse, initial acceleration, etc.). 

The tankage mass of the electric stage m-(- is usually taken to be propor- 
tional to the propellant mass of the electric stage m~. Thus, if k is the 
proportionality constant, 

mt = kmp (2) 

Defining the overall specific powerplant mass a to he the ratio of the 
powerplant mass (including power conditioning and ion engine system) at 1.0 AU 
divided hy the power supplied to the thrustors P and defining the thrustor 
efficiency rj to he the ratio of the propulsive power divided hy the total 
power supplied to the thrustors result in the powerplant mass fraction being 
rewritten as 

^pp p a o I eSo a 

= — a, = - (3) 

m o m o 47.47X10 3 ^ 

where a G is the initial acceleration of the electric stage, I e is the 
specific impulse of the electric stage, and 47.47X10^ is a constant that is 
required for the system of units employed in this paper. 

The payload mass of the chemical stage (initial mass of the electric stage) 
can he written as 

m Q = M 0 - M p - M h 

where the hardware mass is composed of such things as tankage mass. 
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engine mass, guidance and control mass, and structure mass. For chemical 
rockets, the hardware mass can usually he taken to he proportional to the 
propellant mass Mp. If K is the hardware proportionality constant, then 

M h = KM p (4) 

If the chemical stage imparts an impulsive velocity change from circular 
orhit velocity V c to humout velocity the chemical stage payload ratio 

is given hy 


m. 


rf = (1 + K)e 
M n 


-(Vb-VcO/lego 


- K 


(5) 


If equations (2) to (5) are substituted into equation (l) along with the 
final mass relation m^ = m Q - mp, there results 


Mo 


■f 


(1 + K)e 


-(v b -v c )/i c g 0 


-I 


(1 + k) 


“f a o^eSc 


mo 47. 47X1 0 3 


n V m o/_ 


(6) 


This equation is the final form of the function to he maximized. The electric 
stage final mass ratio mf/m D increases as a D , I e , and are increased. 

However, increases in a 0 and I e also increase mpp/m Q , while increasing 
decreases m 0 /M 0 . Clearly, there exist values of a Q , I e , and that 

will result in a maximum payload ratio. These three variables were optimized 
hy a three-dimensional search scheme. 

The maximization of m-p/m 0 can he formulated as a Mayer problem in the 
calculus of variations wherein mission time is treated as a parameter and the 
constraints are composed of the two-dimensional equations of motion plus some 
constraint on the electric engine operation. In the case of constant power 
this latter constraint was chosen to he either maximum thrust or no thrust 
(coasting flight). In the more general case where power is some function G(rj 
of the radius, the constraint was chosen to he either no thrust or thrust 
directly proportional to G(r), where the variation is accomplished hy 
adjusting the propellant flow rate while holding the specific impulse constant. 
(Actually, G(r) is defined to he the ratio of the power at radius r to the 
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power at r = 1. 0 AU. ) The solution to the Euler- Lagrange variational equa- 
tions with an arbitrary power profile G(r) included is given in Appendix A. 
Constant power operation is just a special case of this formulation in which 
G(r) = 1. For this study, the heliocentric travel angle is left unspecified 
(free for optimization). Another boundary condition that is not completely 
specified is the initial heliocentric velocity Vo. which is obtained by 
vectorially adding the velocity relative to the Earth at the sphere of influ- 
ence V s to the Earth's orbital velocity about the Sun V® as shown in 
Fig. 1. The optimum orientation of V s is in the direction of the initial 
electric thrust vector as can be shown by the transversality condition of the 
Mayer formulation. The magnitude of V s is related to the burnout velocity 
V-jj by the conservation of energy as follows: 

v s = yfe - 2vf(l - |2) (7) 

For this study, the sphere of influence radius r s was taken to be 120 Earth 
radii. 

Defining the thrustor efficiency is necessarily somewhat approximate 
because the functional relation Tj(l e ) (e.g. , fig- 2) is dependent on the 
thrustor employed and its state of development. The bulk of the data was 
calculated for electron bombardment thrustors by using the efficiency function 
displayed in Fig. 2. After the constants V c , I e , K, k, and m s /m 0 are 
assigned values (see table I), the payload ratio as given in Eq. (6) may be 
maximized by a three-dimensional search on Vb, a Q , and I e over the class of 
optimal heliocentric trajectories for any given pair of mission time and a. 
Since Vb is constrained to take on values at least as large as the escape 
velocity V e , it might be expected that the optimum Vb could fall exactly on 
this limit (i.e. , Vb jQ pb = V e ) . Ref. 9 shows, however, that the mathematical 
model that has been created for the hybrid system prevents Vb ? opt from 
becoming less than or equal to V e . 
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All electric system. - The payload maximization scheme for the all electric 


system may he performed similarly to the hybrid system scheme. The booster 
burnout velocity Vb is not involved so that only two variables (a D and I e ) 
need be optimized by the search scheme. Optimum Earth- escape trajectories for 
low- thrust vehicles are veiy nearly tangential thrust spirals that can be quite 
accurately computed in closed form. In this case, the ratio m 0 A*o -*- n (-0 

is unity; and the Earth- escape spiral propellant mass and time were computed 
from the equations in Ref. 9 that were developed in Ref. 14. 

The solar cell power function G(r) . - The output power of solar cells 
undergoes changes during the flight due to the following circumstances (l) the 
solar energy flux varies as the inverse square of the radius, (2) the solar 
cell efficiency is a decreasing function of temperature, and (3) the incidence 
angle of the flux with the array surface may change. Normally, the latter 
factor is not particularly inportant, since, for most radii, an array perpen- 
dicular to the flux produces maximum power. But solar cell temperature : 
increases as the radius decreases, and there exists a radius below which the 
temperature function dominates over the flux function, which causes power output 
to fall rapidly. At radii smaller than the maximum power radius, however, the 
array may be tipped at such an angle so as to keep the power constant at its 
maximum value. The tip angle may be increased as radius decreases until the 
back face of the solar panel array is exposed to the Sun (the size of the Sun 
must be considered at small radii) . If the radius is allowed to fall below 
this limiting point, the power will fall rapidly to zero. The power profile 
G(r) used in this study is derived in Appendix B and displayed in Fig. 3. 
Matching the propulsion system load to the changing solar array output power 
and voltage characteristics can be accomplished by changing the number of 
operable thrustors at selected time points and by continual adjustment of the 
propellant flow rate per thrustor. A detailed analysis of this method of power 
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matching is reported in Ref. 13. 

CONSTANT POWER PERFORMANCE 

The payload delivered to 0.1 AU by a Saturn IB/Centaur/electric stage 
vehicle is displayed in Fig. 4 as a function of mission time and specific 
powerplant mass. The family of curves representing the hybrid systems do not 
blend smoothly into the all electric curves because of the restrictive assump- 
tion requiring to be at least escape velocity. (A fuller discussion 

appears in Ref. 9.) The results include three types of electric engine clusters: 
(l) 10- cm- diameter, electron-bombardment engine clusters, (2) 50- cm- diameter, 
electron-bombardment engine clusters, and (3) electrothermal engine clusters. 

The characteristic efficiency curves for these engine types are shown in Fig. 2. 
The gain in efficiency due to increased thrustor size is reflected by a sub- 
stantial performance increase. The performance of the electrothermal engines 
is inferior to that of the electron-bombardment engines because of the unfavor- 
able efficiency curve. Each point on Fig. 4 has been optimized with respect to 
Vfc and I e . Only the circled points include the optimization, with respect to 
a Q . The other points all lie on curves representing all propulsion operation 
(no coast). It is evident that very little payload benefit results from 
coasting trajectories - thus justifying the selection of all propulsion trajec- 
tories throughout much of this paper. Two simplifications arise from this 
constraint: (l) the boundary- value problem associated with the variational 

problem becomes much less sensitive, and, therefore, the iterative method of 
solution converges faster, and (2) a 0 is no longer a variable that needs to be 
optimized} instead, it is determined by the associated boundary- value problem. 

Two typical trajectories of the constant power hybrid system are shown in 
Fig. 5. Arrows denoting optimal thrust direction are placed on the trajectory 
at equal time increments. Trajectory A is of greatest interest since it repre- 
sents an optimal trajectory for the necessarily long mission time trips of high 
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specific poverplant mass systems. It is evident from the trajectory diagram 
that the long missions are characterized by an initial phase that moves the 
vehicle to a high radius (and accompanying low velocity) and a terminal phase 
that simply removes angular momentum. 

SOLAR CELL POWER PERFORMANCE 

The constant power trajectories associated with the systems of major 
interest (a = 50 to 100 lb /kW e ) lie well outside of the Earth’s orbit. Admit- 
tedly, the inclusion of the solar cell power profile presupposes either a major 
performance decrease or a marked alteration of the trajectory to counter the 
effect of a decreasing power profile. The former can be noted in Fig. 6 where 
reduced payloads resulting from using G(r) of Fig. 3 and trajectories differ- 
ing but little from trajectories of constant power are shown by the dashed 
lines. The latter expectation prompted a search for a different class of 
optimal trajectories that would produce higher payloads. Such a class of 
optimal trajectories was indeed discovered. The new class of optimals is 
characterized by inward paths of 2^ revolutions about the Sun. A represent- 
ative trajectory from the 2 ^ revolution class optimals is shown in Fig. 7 where 
thrust pointers have also been included. Varying trip time from 400 to 500 days 
produced only about 0. 2 degree difference in optimal heliocentric travel angle 
(905°). This fact leads to a corollaiy; namely, that a set of optimal classes 
probably exists, each class characterized by the number of revolutions N + 

(N an integer) about the Sun. The class that is globally optimal would be 
primarily a function of trip time. In any event, the present study was restric- 
ted to the 2 tj revolution class optimals. 

Before discussing the payloads associated with these trajectories, a few 
additional comments concerning the nature of these inward trajectories should 
be made. An important parameter in communication studies of solar probes is 
the angle <p between the Sun-Earth radius and the Sun-probe radius. The Sun 
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is a powerful source of background noise, and this creates a communications 
blackout region when cp is near 0° or 180°. The value of cp at perihelion is 
most critical for the solar probe since this is where the major data is to be 
gathered. The values of cp at the three perihelia for the 400 day trip shown 
on Fig. 7 are about 105°, 90°, and 150°, which are perfectly acceptable. 
Increasing the mission time to 450 days would cause all three cp values to be 
within 15° of 90°, an ideal situation from the communications standpoint. 

Another point of interest is the favorable locations of the three optimum 
coast phases. Since these are centered at the perihelia and extend to about 
0. 3 AU, the power normally required by the electric propulsion system becomes 
available to the scientific payload and communications system at a time when 
these systems would benefit most from a substantial power boost. This factor 
can be of considerable importance when comparing electric systems with chemical 
and nuclear systems for missions of high electric power requirements. In such 
cases, electric systems do not need to carry along a separate power supply as 
do the high thrust systems. 

The payloads associated with the 2^ revolution inward trajectories are 
represented as solid curves in Fig. 6. The inward trajectories show a very 
large increase in payload over the outward trajectories. (This is not true in 
the case of constant power. ) A comparison between the optimum constant power 
case (outward trajectories) and the optimum solar cell power case (inward 
trajectories) can be made by using the top four curves in this figure. Clearly, 
neither case is always superior to the other, but as a general rule, solar cell 
power will deliver more payload than constant power at sufficiently high a or 
sufficiently low mission time. In particular, for 400 day missions, solar cell 
power delivers more payload if a is greater than 65 lb/kW g . In addition to 
the Saturn IB/Centaur launch vehicle capability shown in Fig. 6, the payloads 
achievable with an Atlas/Centaur launch vehicle are also shown for the case of 
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solar cell power. The circled point on Fig. 6 represents the coasting trajec- 
tory of Fig. 7 (the curves represent all propulsion) using solar cell power 
and shows that, as in the case of constant power, a small payload "benefit 
results from the use of coast phases. 

According to several recent studies overall solar cell specific 

powerplant masses of 75 lh/kW e are within the realm of current state-of-the-art 
technology. With this in mind, a summary of system performance and design 
parameters for chemical, nuclear, and electric final stages atop Saturn IB/ 
Centaur and Atlas/Centaur launch vehicles is presented in table II for 
a = 75 lb/kW e . The all chemical system data and the 1-i stage nuclear system 
data is presented only at the 75-day mission time since this time represents a 
minimum energy transfer. Neither of these high thrust systems can deliver 
positive payloads. (This data is taken from an unpublished study by the 
Advanced Development and Evaluation Division At the Lewis Research Center and 
Ref. 7.) Likewise, the all electric system (no chemical boost out of the 
initial circular orbit) cannot deliver positive payloads. Only the chemicals 
electric hybrid system can achieve a 0. 1 AU perihelion radius with a positive 
payload. If a larger, more efficient launch vehicle were assumed (e. g. , Saturn 
V) , chemical propulsion alone could deliver some payload^, but by the same 
token, the hybrid system's payloads would increase correspondingly. The 
nuclear rocket system is similar to the chemical system in that larger boosters 
are necessary to provide positive payloads. Besides, nuclear rocket systems do 
not fall into the near future categoiy. Clearly, if both Saturn V class 
boosters and nuclear rockets were available, the nuclear- electric hybrid (instead 
of the chemical- electric hybrid) should be compared with the nuclear system. 

CONCLUDING REMARKS 

To accomplish the 0.1 AU solar flyby mission starting from an initial cir- 
cular Earth orbit, a hybrid system consisting of high-thrust and low-thrust 


14 



stages offers distinct advantages over either system separately. All chemical 
systems are quite unattractive for the close solar probe mission in view of 
their very small payload ratios. All electric propulsion out of a low circular 
Earth orbit is also unattractive unless overall specific powerplant masses of 
less than 50 lb/kW e become available. Electric propulsion systems boosted to 
hyperbolic velocities, however, can deliver significant payloads with the 
relatively high values of a characteristic of the current state-of-the-art 
technology of solar cell power supplies. For instance, a Saturn IB/Centaur/ 
solar electric vehicle with an overall specific powerplant mass of 75 lb/kW e 
could deliver a 1310 lb payload to 0. 1 AU in 400 days, or a 1850 lb payload in 
500 days. 

A more detailed analysis would be required for an in-depth systems design. 
Fig. 7 of the text for exanqple, reveals that a reverse circumferential thrust 
program ( simpler to employ) would be nearly as good as an optimal thrust 
program. The present analysis is sufficient, however, to show that current 
state-of-the-art solar electric propulsion is quite attractive for the 0.1 AU 
solar flyby mission in terms of payload capacity, favorable communications 
angles, and the possible use of the propulsion power supply to act also as 
the power supply of the payload package and communications systems at the 
perihelia. 

APPENDIX A - VARIATIONAL EQUATIONS FOR RADIUS DEPENDENT POWER 

An analysis of the variational problem is not presented here. Instead, 
only the final equations are listed as extensions of those derived in Ref. 15 
for the constant power case. Although the calculations in this paper were 
carried out in two dimensions, the solution of the three dimensional case is 
given here for completeness. The problem is to find the thrust program that 
minimizes the fuel consumption for a transfer trajectory satisfying specific 
initial and final conditions, where the specific impulse is held fixed and the 
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propellant flow rate varies as G(r) if the engines are turned on or is zero if 
the engines are turned off. Let x, y, and z he rectangular position coordi- 
nates and u, v, w the corresponding velocity components in a Cartesian 


coordinate system. Let m he the mass of the vehicle, c the exhaust velocity, 
and Pmax “the propellant flow rate at r = 1. 0 AU for engine- on operation. 

Then the solution to the stated problem is given hy the following set of first 
order differential equations, where differentiation with respect to time is 

denoted hy a superscript dot. 

• c Po / . ^1 

u = - + — G ( r ) X 


m 

c Po 
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r = -\/x^ + 


K = — - Ay 
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c — g 0 I e 

G‘(r) = g 


■•{:: 


if K < 0 


= (W if K > 0 


APPENDIX B - SOLAR CELL POWER PROFILE 


An isolated flat plate solar panel at distance r from the Sun is con- 
sidered where r is large enough to assume that the solar flux lines are 
parallel. If the panel is inclined at an angle i to the solar flux lines 
and is at the equilibrium temperature T, then for equilibrium conditions the 
absorbed power is equal to the sum of the radiated power and the electrical 
output power, that is, 

fe + °f) ^ sln 1 - ( £ c - + -f)^ + " iT) ^ sin 1 (E1) 

where a is the absorptivity, e is the emissivity, z is the packing factor, 
a is the Stefan- Boltzmann constant, I is the solar flux intensity at 1 AU, 
t](T) is the temperature dependent solar cell efficiency, and the subscripts 
c, f, and b refer to cell area, front (toward the Sun) noncell area, and back 
area, respectively. Solving equation (Bl) for T yields 



For temperatures between 150° and 500° K, the last factor of this equation 
ranges between 0.95 and 1.0, and is therefore assumed to be 1.0 in this simple 
analysis. The resulting maximum temperature error is about 10° K at T = 250° K. 
*This derivation was derived by Charles Zola of the NASA Lewis Research Center. 
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The first factor of this equation can "be reduced to a constant hy assigning 
values to the surface properties. For silicon solar cell panels, the following 
values were assumed: Og = 0.94, e c = 0.875, z = 0.95, of = 0. 1, Gf = 0.1, 
et, = 0. 96, I = 130 W/ft 2 , and a = 0.527X10" 8 w/ft 2 /°K 4 . Then equation (B2) 
is simplified to: 

T = 333 (S --' °K (E3) 

Further, Ref. 16 shows that t)(T) is nearly linear over the temperature 
range 150° to 500° K. Thus, if i} 0 is the cell efficiency at T = 298° K, 
n(T) is 

n(T) = n 0 + |g (T - 298) <B4) 

If the data of Fig. 8 in Ref. 16 is used, q 0 = 0. 10, d^/dT = -0. 00046° K _1 , 
and hence, 

tj(T) ® 0.237 - 0.00046 T (B5) 

The electrical output power may now "be calculated hy using Eqs. (B3) and (B5): 


P = n ( T ) -L sin i 


P = 


ii°r o . 


237 - ( S in i) 


V? 


in i) 1 / 4 


sin i 


(B6) 


The equation for the power may he rewritten hy normalizing with respect to 
the maximum power (i = jt/2) at 1.0 AU, P 0 . Thus, if r has units of AU, then 


G(r) 


_ P 2.825 . . 1.825 ^5/4 

= — = — sin 1 7275 ( S3 - n 


(B7) 


(B8) 


Po T c T c 

Solving dG/dr = 0 for sin i yields the optimum value of sin i, 

(sin i) opt = 2. 35 r2 (0 £ r ^ 0. 652) 

(sin i) op t = 3- ( r > 0.652) 

In addition there exists a lower limit on i. At angles less than the lower 

limit, the hack side of the panel is exposed to solar flux (due to the finite 
size of the Sun), which causes the temperature to increase rapidly. The angle 
at which this occurs is given hy tan i = R/r, 


ift 


where R is the radius of the Sun. By combining this expression with Eq. (B8), 


the radius at which this occurs is calculated to be 0. 13 AU. Furthermore, 
fixing i at its lower limit causes G(r) to fall rapidly to zero at slightly 
less than 0.1 AU. Observing these limits on i and substituting Eq. (B8) into 
Eq. (B7) lead to the final form of the solar cell power profile (displayed in 



Fig. 3) : 

* 

(r < 0.13 

( 0.13 £ r £ 0.652) 

(0.652 £ r) 

-£**-< I**-* * 

It should "be realized that this derivation is somewhat approximate for a number 
of reasons. For instance, experiments performed on silicon solar cells-*- ^ 
indicate that cell efficiency depends on solar flux intensity and that "both 
cell efficiency and absorptivity are also functions of the incidence angle i 
when i is less than about #0°. Also, performance degradation due to charged- 
particle irradiation-*-® and the interaction between the spacecraft and the solar 
panels were also ignored. These effects are not as yet defined well enough to 
warrent inclusion in this preliminary analysis. The silicon solar cell power 
profiles appearing in Refs. 17' and 18 do not differ greatly from Fig. 3. In 
Ref. 18 gallium- arsenide solar cell power profiles were found to be considerably 
more favorable than those of silicon cells for close-in solar probes. 
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TABLE I. - ASSUMED CONSTANTS 


Parameter 

Assumed value 

V c , ft /sec 

25 600 

mg/nio 

0. 10 

k 

0.10 (electron-bombardment engines) 
0.15 (electrothermal engines) 

K a 

0. 137 ( Saturn IB/Centaur) 
0.447 (Atlas/Centaur) 

I c , sec a 

420 (Saturn IB/Centaur) 
440 (Atlas/Centaur) 


a The values of K and I c appearing in this tahle 
should not he taken as numbers corresponding to 
actual Centaur hardware. Instead, they are two 
parameters evaluated for a curve fit of the form 
given by Eq. (5) to Saturn IB/Centaur and Atlas/ 
Centaur launch vehicle performance data. 


TABLE H. - SYSTEM PERFORMANCE FOR 0. 1 AU SOLAR FLYBY MISSION 8, 


System 

Mission 

time, 

-k* 

days 

Payload, 

lb 

Power to 
thrustors, 
k¥ e 

Specific 
impulse 
of final 
stage, 
sec 

High 

thrust 

burnout 

velocity, 

ft/sec 

Saturn IB/Centaur/kick stage 

75 

0 

— 

444 

67 000 

Saturn IB/nuclear (l-| stages) 

75 

0 

— 

800 

67 000 

Saturn IB/Centaur/electric 

400 

1050 

48.0 

3900 

39 200 

(constant power) 

500 

2350 

49.9 

4400 

37 800 

Saturn IB/Centaur/electric 

400 

1310 

31.9 

4100 

41 800 

( solar cell power) 

500 

1850 

37.1 

4900 

39 900 

Atlas/Centaur/electric 
(solar cell power) 

500 

230 

9.4 

4700 

37 600 


Saturn IB assumed to inject 32 000 lb payload into 100-n. mi. Earth orbit. 
Atlas/Centaur assumed to inject 10 800 lb payload into 100-n. mi. Earth orbit. 
Overall specific powerplant mass a = 75 lb/kW e . 
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Fig. 5. - Two typical 0. 1 AU solar fly-try trajectory diagrams for constant power 
hybrid system. Optimum V^, 8, 


Assumed initial mass in 
100 n. mi. Earth orbit 



Fig. 6. - Payload delivered to 0. 1 AU by Saturn IB/Centaur/electric 
vehicle using 50 cm electron bombardment engines. 
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